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LONG-TERM  GOAL 

My  long-term  goal  is  to  make  substantial  contributions  to  enhancement  of  forecast  skills  concerning 
the  oceans  and  atmosphere,  by  deepening  our  knowledge  for  the  nature  of  predictability.  I  place  my 
emphasis  on  transition  mechanism  between  dynamical  regimes  of  planetary  flows  and  coherent 
structures  from  both  Eulerian  and  Lagrangian  points  of  view. 

OBJECTIVES 

I  wish  to  develop  theoretical  frameworks  for  enhancement  of  predictability  based  on  dynamical 
systems  theory.  To  help  enhance  predictability  of  sudden  and  severe  events,  I  wish  to  investigate 
transition  mechanisms  between  dynamical  regimes.  I  also  wish  to  throw  a  bridge  between  Eulerian 
and  Lagrangian  viewpoints  and  explore  impact  of  severe  events  on  large-scale  transport  in  planetary 
flows.  Knowledge  obtained  by  these  predictability  studies  will  lead  to  a  design  of  comprehensive 
ensemble-forecast  systems  and  data-adaptive  observing  systems. 

APPROACH 

To  achieve  my  goals,  I  plan  to  develop  new  theories  and  improve  already  existing  methodologies  so 
that  they  can  be  combined  systematically.  I  start  from  assessing  and  extending  individual  elements  of 
the  theoretical  framework,  and  develop  new  theories  to  fill  the  gaps  between  them  as  necessary.  My 
approach  involves:  identification  and  detection  of  the  predictability  elements,  Eulerian  and  Lagrangian 
descriptions  of  the  probability  density  function  evolution,  treatment  of  nonlinearity  as  well  as  additive 
and  multiplicative  stochasity  in  data  assimilation  systems.  To  improve  forecast  skills  for  severe  and 
sudden  events,  I  investigate  role  played  by  stochastic  noises.  Severe  events  can  be  viewed  as  rare 
extreme  bursts  and  therefore  may  be  related  to  stochastic  noises  whose  probability  distribution  has 
heavy  tail.  I  construct  a  new  innovative  methodology  for  data  assimilation  systems  subject  to 
dynamical  and  observational  noises  with  heavy-tail  distributions.  To  explore  impact  of  transient  flow 
dynamics  on  large-scale  transport,  I  combine  geometrical  approach  of  dynamical  systems  theory  and 
spatio-temporal  data  analysis  technique. 
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WORK  COMPLETED 


An  innovative  method  for  direct  assimilation  of  Lagrangian  data  in  meteorology  and  oceanography  has 
been  formulated  based  on  the  extended  Kalman  filter  (EKF).  The  method  has  been  applied  to  the  point 
vortex  systems.  A  classification  of  the  Lagrangian  trajectories  has  been  made  for  practical  applications 
of  the  direct  Lagrangian  data  assimilation  method  to  complex  ocean  circulation  models. 

In  order  to  examine  the  predictability  of  the  coherent  structures  in  the  stratified  ocean,  a  layer-model 
for  point  vortex  system  has  been  investigated.  Judicious  observing  system,  including  parameter- 
estimation  capacity,  has  been  designed  for  improved  estimation  of  the  flow  dynamics  in  the 
unobserved  sub-surface  layers. 

To  study  predictability  concerning  transition  from  oscillatory  to  singular  behavior,  an  idealized,  highly 
nonlinear  dynamical  system  has  been  developed.  The  model  is  based  on  interplay  between  restoring 
force  and  hysteria  in  inertia,  which  can  be  found  in  many  physical  systems.  The  bifurcation  diagram 
has  been  investigated  to  enhance  forecast  skill  for  the  sudden  transition. 

A  hybrid  transport  theory  (TIME:  transport  induced  by  mean-eddy  interaction)  has  been  formulated  to 
estimate  transport  induced  by  the  anomaly  in  the  large-scale  atmosphere  and  ocean  flows.  It  has  been 
extended  to  identify  signature  and  analyze  role  of  variability  in  application  to  a  large-scale  atmospheric 
flow. 

A  comprehensive  study  to  analyze  transport  in  the  Monterey  bay  surface  flow  has  been  pursued.  In 
contrast  to  conventional  Eulerian  analysis  techniques,  Lagrangian  synoptic  maps  have  been  used  to 
study  the  Lagrangian  ocean  surface  properties. 

To  address  extended-range  predictability  in  terms  of  the  multiple-regime  paradigm  and  preferred 
transition  paths  between  them,  the  three-level  quasi-geostrophic  (QG)  model  has  been  investigated. 
Preferred  transitions  paths  are  examined  also  for  a  stochastically  forced  Lorenz  system,  to  help  explain 
the  striking  feature  of  the  QG  model. 

RESULTS 

The  new  Lagrangian  data  assimilation  system  augments  the  state  variables  with  the  Lagrangian  tracer 
variables  in  addition  to  the  prognostic  variables  of  the  dynamic  model.  The  EKF  computes  not  only 
the  state  variables  but  also  the  error  covariance  matrix.  Through  the  cross  correlation,  the  new 
assimilation  system  updates  the  unobserved  prognostic  variables  using  the  observed  Lagrangian 
tracers.  This  is  based  on  the  same  idea  as  the  parameter  estimation  using  the  EKF.  Figure  1  shows  a 
successful  application  of  the  EKF  to  the  surface  layer  observing  system  for  estimation  of  the  dynamics 
using  the  point-vortex  model  that  is  the  lowest-order  representation  of  the  coherent  structures  in  the 
large-scale  ocean  flows.  Performance  of  the  data  assimilation  system  crucially  depends  on  the 
judicious  choice  of  the  key  parameters  in  the  observing  system,  in  particular  observation  frequency. 
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Figure  1.  Surface  layer  observation  system  for  subsurface  layer  estimation. 

[Left:  underlying  truth  by  control  run  with  surface  (blue)  and  sub-surface  (red)  ocean  dipole;  Right: 

successful  data  assimilation  with  parameter  estimation] 


By  applying  the  mixture  model  clustering  technique  to  leading-order  empirical  orthogonal  functions  of 
the  500hPa  level  streamfunction  field  obtained  by  the  QG  model  simulation,  the  four  robust  regimes 
roughly  corresponding  to  opposite  phases  of  the  Arctic  Oscillation  and  the  North  Atlantic  Oscillation 
are  found  (Figure  2).  The  preferred  transition  paths  are  found  not  to  be  aligned  with  the  line  segments 
between  regime  centroids  in  most  cases.  It  might  point  to  heteroclinic  or  homoclinic  connections 
between  unstable  equilibria  as  demonstrated  using  the  Lorenz  equations. 

IMPACT/APPLICATION 

The  new  data  assimilation  methodology  currently  being  developed  in  this  work  can  fill  the  gap 
between  data  assimilation  and  dynamical  systems  approach  to  Lagrangian  transport  theory,  which  has 
been  one  of  the  focuses  in  the  ONR  predictability  DRI  program.  The  methodology  is  general  and  can 
be  applied  to  numerous  observations  that  have  not  been  assimilated  using  the  conventional  assimilation 
methods. 

Extended-range  weather  prediction  depends  in  a  crucial  way  on  skill  at  forecasting  the  duration  of  a 
regime  event  or  other  persistent  anomaly,  as  well  as  the  subsequent  onset  of  another  persistent  anomaly 
after  the  break  of  the  current  one.  The  result  obtained  may  suggest  that  a  similar  analysis  of  numerical 
weather  prediction  models  could  help  forecast  well  in  advance  regime  breaks  and  subsequent  onsets. 

To  do  so  would  require  an  efficiently  designed  observational  system  that  is  able  to  track  spatial 
signatures  of  regime  transitions. 

Extending  TIME  theory  offers  a  variety  of  new  directions  for  large-scale  planetary  flows,  including 
analysis  of  transport  processes  and  mechanism,  impact  of  specific  events  in  the  flow  dynamics,  as  well 
as  role  played  by  variability. 
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Figure  2.  Streamfucntion  filed  of  the  four  atmospheric  weather  regimes  at  500hPa. 
[Positive  (+)  and  negative  (-)  phases  of  North  Atlantic  Oscillation  (NAO)  and  Arctic 

Oscillation  (AO)] 
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